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ABSTRACT 



We present quantitative investigations of the weak lensing effect on the two-point correlation 
functions of local maxima (hotspots), ^pk-pk(^), in the cosmic microwave background (CMB) 
maps. The lensing effect depends on the projected mass fluctuations between today and the 
redshift Zrec ~ 1100. If adopting the Gaussian assumption for the primordial temperature 
fluctuations field, the peak statistics can provide an additional information about the intrinsic 
distribution of hotspots that those pairs have some characteristic separation angles. The weak 
lensing then redistributes hotspots in the observed CMB maps from the intrinsic distribution 
and consequently imprints non-Gaussian signatures onto ^pk-pk(^)- Especially, since the 
intrinsic ^pk-pk(^) has a pronounced depression feature around the angular scale of ~ 70' 
for a flat universe, the weak lensing induces a large smoothing at the scale. We show that the 
lensing signature therefore has an advantage to effectively probe mass fluctuations with large 
wavelength modes around A « 50/i~^Mpc. To reveal the detectability, we performed numerical 
experiments with speciflcations of MAP and Planck Surveyor including the instrumental effects 
of beam smoothing and detector noise. It is then found that our method can successfully provide 
constraints on amplitude of the mass fluctuations and cosmological parameters in a flat universe 
with and without cosmological constant, provided that we use maps with 65% sky coverage 
expected from Planck. 



Subject headings: cosmology: theory - cosmic microwave background - gravitational lensing - 
large-scale structure of universe 



Introduction 



The temperature anisotropies in the comic microwave background (CMB) contain detailed information 

The recent high precision 



about the underlying cosmological model ( Hu, Sugiyama fc Silk 1997 ) 
balloon-borne experiments. Boomerang ( pe Bernardis et al. 2000| ; |Lange et al. 2000D and MAXIMA 



(Banany et al. 2000; Balbi ct al. 2000), revealed that the measured angular power spectrum C; are in 
good agreement with that predicted by standard inflation paradigm. On the other hand, the large-scale 
structure of the universe imprints secondary effects on the primordial temperature fluctuations. One of 
them is the weak lensing effect; the CMB photons are randomly deflected by the gravitational field due to 
the intervening large-scale structure during the propagations from the last scattering surface to a telescope. 
The weak lensing can be a powerful probe for mapping inhomogeneous distribution of dark matter in the 



universe (Gunn 1967; see Bartelmann & Schneider 2000 for a review), which is not directly attainable 



by any other means. In fact, several groups (Van Waerbeke et al. 2000; Bacon, Rcfregier & Ellis 2000 



Wittman et al. 2000; Kaiser, Wilson & Luppino 2000) recently have reported significant detections of 
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the coherent distortion of faint galaxies images arising from the gravitational lensing by the foreground 
large-scale structure, and showed that those results can provide some constraints on the cosmological 
parameters. However, it would be still extremely interesting to be able to measure the lensing effects on the 
CMB and the detection would be very precise for constraining the cosmological parameters because there 
is no ambiguity in theoretical understanding of the unlensed CMB physics and about the distance of the 
source plane. 

The inflationary scenarios also predict that the primordial temperature fluctuations are Gaussian. 
In this case, the statistical properties of any unlensed CMB field can be exactly predicted based on the 
Gaussian random theory developed by Bardeen et al. (1986; hereafter BBKS) and Bond & Efstathiou (1987; 
hereafter BE) for three- and two-dimensional cases, respectively. However, the weak lensing then induces 
the non-Gaussian signatures in the observed CMB maps. Based on these considerations, some specific 
features on the lensed temperature fluctuations field have been revealed. Bernardeau (1998) investigated 
how the weak lensing alters the probability density function (PDF) of the ellipticities defined from the local 
curvature matrix of the temperature fluctuations field. The unlensed PDF indeed has specific statistical 
properties for the two dimensional Gaussian field, and then the gravitational distortion induces an excess of 
elongated structures in the CMB maps in the similar way as the distortions of distant galaxies. Although 
the method could be a powerful probe to measure the lensing signatures around the characteristic curvature 
scale (~ 5') of the unlensed temperature field, the instrumental effect of a finite beam size is crucial for 
the detection because the beam smearing effect again tends to circularize the deformed local structures. 
Hence, Van Wearbeke, Bernardeau & Benabed (1999) investigated how the weak lensing causes a coherent 
distribution of the relative orientation between the CMB and distant galactic ellipticities, and proposed 
that it can be a efficient tool because the orientation of the CMB ellipticities is robust against the beam 
smearing. 



We ( Takada, Komatsu fc Futamase 2000 ; hereafter TKF) recently investigated the weak lensing effect 
on the two-point correlation function of local maxima (hotspots), say ^pk-pk(^), in the two dimensional 
CMB maps. Since the distribution of hotspots is a point process, the analysis focused on the secondary 
effect how the weak lensing redistributes hotspots in the observed CMB maps from the intrinsic distribution. 
The unlensed ^pk-pk(^) can be then accurately predicted by the Gaussian random theory once Ci is given 



(BE; Heavens & Sheth 1999 hereafter HS). According to the acoustic peaks in the C;, the pairs of hotspots 
have some characteristic angular separations on the last scattering surface. We then found that the weak 
lensing fairly smooths out the oscillatory shape of ^pk-pk(^)- In particular, the most interesting result is 
that the lensing contribution to ^pk-pk(^) at angular scales (w 70') corresponding to the first Doppler peak 
of Ci is relatively large, and we thus expect that Cpk-pk(^) can be a sensitive statistical tool to measure 
the projected mass fiuctuations at such larger scales than the other methods do. The crucial quantity 
of our method is the lensing fiuctuations of relative angular separation between two CMB photons, and 
the lensing signatures to ^pk-pk(^) at such large scales is the consequence of large scale modes of lensing 
deflection angles. The simulated maps indeed illustrate that each displacement of peak positions in the 
lensed maps from the unlensed maps is relatively large even though both global features of pattern of 
temperature fluctuations nearly trace each other (Zaldarriaga 200C and also see Figure Furthermore, 



these considerations lead to the expectation that our method is not particularly affected by the beam 
smearing effect. 

The purpose of this paper, therefore, is to quantitatively investigate the weak lensing effect on 
Cpk-pk(^) and reveal in detail the physical interpretations of the effect. For the practical purpose we 
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perform quantitative investigations of the problem whether the future sateUite missions, MyUQ and Planck 
Surveyo"^, can measure the lensing signatures to ^pk-pk(^) for constraining the cosmological parameters. 
This can be done by using numerical experiments of both unlensed and lensed CMB maps including the 
instrumental effects of finite beam size and detector noise. 

This paper is organized as follows. In the next section, we introduce the formalism to investigate the 
lensing effect on ^pk-pk(^^)- In Section 4, the formalism is applied to some specific cosmological models, 
and we then compute the signal-to-noise ratios of the lensing signatures to ^pk-pk(^) using the numerical 
experiments in Section 5. Finally, in Section 6 we present discussions and conclusions. 

2. Weak Lensing EflFect on Two-point Correlations of Hotspots in the CMB: Formalism 

A bundle of CMB photons is randomly deflected by the inhomogeneous matter distributions of the 
intervening large-scale structure as it propagates from the last scattering surface to a telescope. The two 
CMB bundles observed with a certain angular separation thus have a different angular separation when 
emitted from the last scattering surface. The ensemble averages of the second moment of the relative 
separation fluctuations and the following characteristic function can be then easily calculated by using the 
power spectrum approach developed by Seljak (1996); 

(exp[zZ • {50, - 6e2)])^e,-e,\=e - 1 - ^ [4l,oW + cos(2(p04l,2(^)] , (1) 

where S9i{= 66{9i)) and (502(= 50{02)) are the angular excursions of the two bundles, and ( )g 
observationally means the average performed over all pairs with a fixed observed angular separation 
^- fGL.ol^) and (7gl,2(^) represent isotropic and anisotropic contributions to the lensing dispersion, 
respectively. Although the anisotropic one is ignored in TKF for simplicity, we also take into account the 
contribution in this paper. It is convenient to express those dispersions in terms of the logarithmic angular 
power spectrum of the deflection angle, Pgl(Oj ^ 

4l,2(^) = ^ JjPGL{l)J2iW), (2) 

with 

Pgl(0 = 9H^^l^ dxa-^{T)W\x, X.oc)P6 [k = ■ (3) 

The statistical properties of the lensing effects are thus entirely determined by Pgl(0; because the lensing 
field is expected to be also Gaussian at relevant angular scales {9 ^ 10'). r is a conformal time, x = tq — t, 
Jn{x) is the Bessel function of order n, and the subscript and "rec" denote values at present and a 
recombination time, respectively. Ps{k,T) is the power spectrum of matter fluctuations field, Hq = lOO/i 
km s~i Mpc^i and flmo denote the present Hubble parameter and the present energy density of matter, 
respectively. r(x) is the corresponding comoving angular diameter distance, defined as K^^^^ sin K^^^Xj 



^ http : //map . gsf c . nasa . gov 

^ http : //astro . estec . esa.nl/SA-general/Projects/Planck/ 
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X, {—K)~^/'^s\-a\\{—KY/'^x for K > Q^K = Q,K < Q, respectively, where the curvature parameter K is 
represented as K = (Omo + ^ao ~ ^)H^ and ri^o is the present vacuum energy density. The projection 
operator W{x,Xrcc) on the celestial sphere is given by W(x,Xrcc) = f{Xicc — x)/^(Xroc)- In the derivation 
of equation (|^), we have employed two approximations. First is the flat-sky approximation where the two 
dimension Fourier transformation is used neglecting the curvature of the celestial sphere. This is based 
on the consideration that the lensing are important only on small angular scales. Second is the Born 
approximation that the integral can be evaluated along the unperturbed nuU-geodesics of CMB photon. 
Hu (2000) recently investigated the correction to the flat sky approach by directly evaluating the lensing 
effect in harmonic space. This correction to our method, however, is small as will be discussed later. The 
Born approximation is valid as long as the (tqi^(9)/9 <C 1 is satisfied, and we numerically confirmed this on 
the relevant angular scales for all cosmological models considered in this paper. Importantly, magnitude of 
the lensing dispersion (|l|) is particularly sensitive to and the normalization of matter power spectrum, 
which is conventionally expressed in terms of the rms mass fluctuations of a sphere of 8/i^^Mpc, i.e., erg. 

The hotspots are local maxima in the two dimensional CMB sky map, and hence the distribution 
obeys a point process. Once the angular power spectrum of the temperature fluctuations field Ci is given, 
the Gaussian assumption allows us to exactly predict statistical properties of the intrinsic distribution of 
hotspots following the methods developed by BBKS and BE. Since C; has oscillatory features such as series 
of Doppler peaks, the pairs of hotspots are distributed with some characteristic separation angles on the last 
scattering surface as shown by HS and TKF. This result leads to the following expectation. Let us consider 
all pairs of hotspots separated with the certain characteristic angular scale. Although all those pairs 
should be observed with the characteristic scale in the absence of the lensing, they are actually observed 
with various different separations in random lines of sight because of the weak lensing. The probability 
distribution of observed separation angles then has the lensing dispersion at the characteristic scale. 
Since the effect can be measured in only a statistical sense, in this paper we focus on investigations of the 
lensing effect on the two-point correlation function of hotspots. Note that our method do not consider 
spurious hotspots created by the lensing, but it is a good approximation because it has been shown that 
an additional power of the anisotropics generated by the weak lensing is very small and important only at 



small scales (6* < 1') ( |A/Ietcalf fc Silk 1997| ; [Seljak fc Zaldarriaga 2m\ , [Zaldarriaga 200G| ). These features are 
indeed illustrated by the numerically simulated CMB maps (see Figure]^) 

To calculate the weak lensing effect on the two-point correlation function of hotspots in the CMB 
maps, we first define the number density fluctuations field of hotspots as 

Sr^Me) = rfl^, (4) 

where rip]^{9) and npk are the number density field and the mean number density of hotspots above a 
certain threshold v, respectively. The threshold is conventionally expressed in units of the rms temperature 
fluctuations as = Apk/tro, where Apk is the value of temperature fluctuation at the peak position defined 
by A(0pk) = ST{6p]^)/TcMB and the dispersion is defined by CTq = {A^{6)). Similarly, the other spectral 
parameters are defined by = ((VA)^) and (t| = ((V^A)^). These parameters can be expressed in terms 
of Ci in the context of the small angular approximation (BE) as = J {Idl / {2tt))CiP" . The analytical 
expression of npk(> ly) has been derived by BE and is in detail presented in appendix Because of the 
mapping effect due to weak lensing, the lensed (observed) fluctuations field, Sn^^{9), at a certain angular 
position 9 is the intrinsic field at another position 9 + 59 on the last scattering surface, where 50 is the 



- 5 - 



deflection angle. Thus Snp^{9) can be expressed as 



(5) 



where Sn^ is the Fourier component of unlensed field Sn{6). Therefore, since the lensing deflection angle 
induced by the large-scale structure and the CMB field on the last scattering surface are statistically 
independent, the lensed (observed) two-point correlation function of hotspots, C^k^pk(^)i can be calculated 
with help of equation (]^) as 



(27r)2 J (27r)2 

ZTT 



{5ni5ny)U 



i-^4L,o(^)Moae) 



i{l-s9i-l'-502) 



where we have used the following Gaussian random property of 6nj^ ; 



(6) 



(7) 



-Ppk-pk(0 is the angular power spectrum of unlensed two-point correlation function of hotspots, ^pk-pk(^)- 
The derivation of ^pk-pkl^*) is presented in detail in the appendix^. The relation between ^pk-pk(6') and 
^pk-pk(0 can be expressed as 

Fpk-pk(0 - 27r / ^?d0epk-pk(0)Jo(^e). 



By using this equation, equation 



can be therefore rewritten as 



1 



(8) 



(9) 



This is the equation which we use for theoretical predictions of the lensing effect on ^pk-pk(^)- If we ignore 
the anisotropic lensing dispersion fTGL,2, the expression (^) can be further simplified (see appendix 
Importantly, this equation indicates that the lensing contribution to ^pk-pk(^') at a certain scale 9 arises 
only from the lensing dispersion ctgl [0] at the same scale. Hence, detections of scale dependences of the 
lensing signatures to ^pk-pk(^) allow us to reconstruct the lensing dispersion at the respective scales, more 
interestingly, to reconstruct the projected matter power spectrum. 



3. Cosmological models 



To make some quantitative predictions, one needs to specify cosmological models. For this reason, we 
adopt following adiabatic cold dark matter models with fimo = Ij ^^Ao = 0, /i = 0.5 (hereafter SCDM) and 
f^mo = 0.3, VL\q = 0.7, h ~ 0.7 (hereafter LCDM), respectively. These models are motivated by the fact 
that the recent high precision measurements of Ci supported a flat universe under the adiabatic condition 



as suggested by standard inflationary scenarios ( de Bernardis et al. 200C ; Hanany et al. 2000 ). The baryon 
density is chosen to satisfy f2bo^^ — 0.019, which is consistent with values obtained from the measurements 
of the primeval deuterium abundance ( Buries fc Tytler 199§| ). As for the matter power spectrum, we 
employed the Harrison-Zel'dovich spectrum and the BBKS transfer function with the shape parameter from 
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Fig. 1. — The dependence of ctgl [d) /O on 9 is shown for two cosmological models with various normalizations 
of CDM power spectrum and ilbo^^ = 0.019. The solid lines show results of erg = 0.5, 1.0 and 1.5 in SCDM 
(fimo = 1 and h = 0.5) from bottom to top while the dashed lines are erg = 1.0,1.5 and 2.0 in LCDM 
(f2mO = 0.3, f2Ao — 0.7), respectively, where the nonlinear effect on the matter power spectrum is computed 
using the fitting formula by Peacock & Dodds (1996). The thin dot-dashed line is for the linear power 
spectrum in SCDM with = 1.5. 



Sugiyama (1995). The free parameter in each model is only the normalization of the present-day matter 
power spectrum, i.e., trg- The nonlinear evolution of the power spectrum can be modeled using the fitting 
formula given by Peacock & Dodds (1996). To compute the angular power spectrum C;, we used helpful 
CMBFAST code developed by Seljak & Zaldarriaga (1996). 



4. Theoretical Results 

4.1. Dependence of lensing dispersion on cosmological parameters 

In Figure |l| we plot the relative lensing dispersion (7Gh{0)/d as a function of the separation angle 6 for 
different sets of erg in SCDM and LCDM models, which are computed using equations (|l]). The dependence 
of (Tg in each model is demonstrated by choices of cg = 0.5, 1, and 1.5 (solid lines) in SCDM, and erg = 1, 1.5 
and 2.0 (dashed lines) in LCDM from bottom to top, respectively. The normalizations from the COBE 
4-year measurements (e.g. Bunn &: White 1997 ) and the X-ray cluster abundance ( Eke, Cole fc Frenk 199fi| ; 
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Fig. 2. — The logarithmic angular power spectrum of the Icnsing dispersion as a function of I for SCDM 
(solid) and LCDM (dashed) models with erg = 1.5 as in Figure |l|. The power spectrum is computed using 
equation (|^). The thin and thick lines in each model show the corresponding results of using linear and 
nonlinear three dimensional power spectrum of matter fluctuations, respectively. 



Kitayama fc Suto 19971) roughly correspond to as = 1.2 and 0.5 for SCDM, and cts - 1 - 1-5 for LCDM, 



respectively. Furthermore, the recent several measurements of cosmic shear (e.g. Van Waerbeke et al. 200C ) 
have suggested as ~ 1.5 ± 0.5 for the current favored LCDM models, while uncertainties involved in redshift 
distributions of distant galaxies, the cosmic variance of the variance of the shear and the systematic errors 
of the signals still remain unresolved. Figure |l| clearly shows that the magnitude of (Tgl(^) at a certain 
angle has a strong dependence on the amplitude of as in each cosmological model. Since the logarithmic 
angular power spectrum (|^) of ctgl has contributions from the proportional factor JI^q and the distance 
and the growth factors, the combination yields a dependence of fimo on ctgl(^). Although the Hubble 
parameter h affects ctgl mainly through the shape parameter F « ^Imoh in the matter power spectrum, the 
dependence is weaker. The thin dot-dashed line shows a result of using the linear matter power spectrum 
for as = 1.5 in SCDM, and it reveals that the effect of the nonlinear evolution on ctgl is not important at 
9 ^ 10', where the two-point correlations function of hotspots has most power of correlations. It is therefore 
expected that the measurements of ctgl can generally provide constraints on as — i^mo plane. To break the 
degeneracy between ft^Q and CTg, we have to measure the scale dependence of ctgl(^') with respect to 9. 

As shown by equations ^ and (^), the Icnsing dispersion arises from the projected matter power 
spectrum Ps{k) weig 

hted with a"2^2^^^^^^^)_ pj 

gure H shows the logarithmic angular power spectrum 
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Fig. 3. — (a) The logarithmic contribution to the 2D power spectrum Pgl(0 defined by equation (^) as a 
function of fc for I — 10, 100, 1000, where the normahzation is arbitrarily scaled. The models are same as 
those in Figure 0. (b) The logarithmic contribution to Pgl(0 for each I mode as a function of 1 + z. 



^gl(0 defined by equation (|^) as a function of /. The thin and thick curves are results of using the linear 
and nonlinear matter power spectra, respectively, and the figure demonstrates that the nonlinear effect is 
important on I }t 1000, which corresponds to angular scales of 6* ^ 20' from the relation of I w 2tt/9. The 
shape of Pql peaks around / « 100 and thus the contributions to the deflection angle of each CMB photon 
come mainly from modes with such large / (see Figure The essential quantity of lensing effect on ^pk-pk 
is lensing fluctuations of relative separation angle, and the term including Bessel function in <7gl,o(^) of 
equation (||) indicates that the lensing fluctuations for two CMB bundles separated with a certain angle 9 
arises dominantly from the integrations of -Pgl(0 over I > 2ti/9 in /—space. This physically means that 
the gravitational lensing effect on the relative angular separation is caused mainly by the projected matter 
fluctuations lied between the two bundles. Accordingly the two bundles with smaller separation are more 
strongly affected by the smaller scale structures of the universe. Furthermore, from these interpretations we 
can conclude that corrections to the flat sky approximation proposed by Hu (2000) do not affect our results 
at / ^ 100 because the corrections are important only at I < 10. 

Since the more fundamental quantity is the three dimensional density fluctuations characterized by 
the power spectrum Ps{k), we have to see the relation between the two power spectra, namely Ps{k) and 
Pgl(0- Figure ^(a) shows the logarithmic contribution to Pgl(0 foi' ^ given Z-mode as a function of three 
dimensional wavenumber k in the same models as in Figure |[ where normalizations of those curves are 
arbitrarily scaled. These functions have relatively broad shape and peaks at A = fc/27r « 630/i^^Mpc for 
Z = 10, at fc « 125/i-iMpc for I = 100, and k « 21/i-iMpc for I = 1000 in the considered SCDM model, 
respectively. These relations between I and k depend on the shape of matter power spectrum. On the other 
hand, since in the LCDM model the comoving angular diameter distance at a certain redshift is larger than 
the corresponding distance in SCDM, the peak-wavelengths are smaller than those in the SCDM cases for 
each I mode. The figure also demonstrates that the higher I mode is affected more strongly by the matter 
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fluctuations with smaller wavelengths. 

The next question is the redshift distribution of the contribution to a given Pgl(0- Figure |^(b) shows 
the logarithmic contribution to -Pgl(0 foi' each I mode as a function of (1 + z). Since the window function 
W/a in Pgl(0 ^ monotonous decreasing function with respect to z and has no characteristic redshift, the 
question which redshift structures give the dominant contribution depends on the shape of matter power 
spectrum. The figure clearly demonstrates that each curve peaks at a certain redshift; for low / mode the 
contribution is dominated by the low-z structures and, on the other hand, high / mode has wide range 
contributions in the redshift space. These results can be explained as follows. As shown in Figure |^(a), there 
is the peak scale fc,nax of three dimensional mass fiuctuations that provides most dominant contribution to 
^gl(0 fo'' each I mode. Because of the projection effect on the celestial sphere, the contribution from the 
fcmax mode comes from structures at a specific redshift Zmax which satisfies the relation of fcmax ~ l/f{zms^x)■ 
As a result, for lower I mode the contribution will be dominated by lower z structures. On the other hand, 
for / = 1000 mode the contribution peaks at z « 0.5 with a long tail expanding to higher z, where the tail is 
caused by the fact that the window function W/a has larger power at lower z. In fact. Figure §(b) confirms 
these interpretations. We could therefore directly probe dark matter clustering at low and high rcdshifts for 
low and high I modes, respectively, in principle by using the weak lensing signatures. 



4.2. Unlensed two-point correlation function of hotspots with the instrumental effects 

The unlensed two-point correlation functions of local maxima (hotspots), say ^pk-pkl^*), in the CMB 
maps can be accurately predicted based on the Gaussian random theory. The derivation is presented in 
the appendix ^ in detail. To show the shape of ^pk-pk(^), we performed a six dimensional numerical 



integration using equation ( p35| ) because the two of the eight dimensional integration in equation (B5) 
can be analytically done. For the practical purpose, we also take into account instrumental effects of 
finite beam size and detector noise. The beam smearing effect on the temperature fluctuations fleld can 
be modeled by the Gaussian beam approximation characterized by a filter function Fi = cxp[—P9'^/2] in 
Z-space, where the smoothing angle Og is expressed in terms of the full-width at half-maximum angle 0fwhm 
of a telescope as 0s = 6'fwhm/V81n2. The noise level of detectors is conventionally expressed in terms of the 
temperature fluctuations per a pixel on a side of FWHM extent as Anoise = cTscns- If we assume that the 
primordial temperature and the noise fields are uncorrelated, by modifying the angular power spectrum 
Ci to Ci = [Ci + cTsons^fwhm) ^'^pI"^^^^]' thcse instrumental effects can be approximately included into the 
theoretical predictions (HS), because in the Gaussian random theory C; contains complete information 
about statistical properties of any intrinsic CMB field. The numerical experiments indeed show that this 
treatment works well. 

Figure ^ shows the unlensed ^pk-pk(^) with and without the experimental effects as a function of 
separation angle 9 in the SCDM and LCDM models, where we considered the two hotspots of height 
above the threshold v — l(Apk = co)- As for the instrumental effects, we have employed the specifications 
of Planck 217GHz channel and MAP 90GHz channel; the beam size and noise level are assumed to be 
Of^Yim = 5.5' and (iscns = 4.3 x 10^^ for Planck while flfwhm = 12.6' and o-gons = 2.5 x 10^^ for MAP. 
The solid lines in each panel show ^pk-pk(^) including both effects of the beam size and noise, and the 
dotted lines in the upper left and bottom left panels show an ideal case without those effects. The ideal 
cases demonstrate that the intrinsic ^pk-pk(^) has a prominent peak at « 10' and the damping tail at 
6 < 10'. These features physically mean that the primary temperature field has the characteristic curvature 
scale of the order of 10', which can be estimated by f?* = \f2a\la2 ^ 5' in Appendix and has smooth 
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Fig. 4. — The unlensed two-point correlation function of hotspots of height above threshold ly — 1, 
Cpk-pk(^), as a function of the separation angle 6 for SCDM (upper panels) and LCDM (bottom panels). 
The instrumental effects of finite beam size and detector noise for the expected Planck and MAP surveys 
are taken into account. The solid and dot-dashed lines shows the cases including both those effects and only 
the beam smearing effect, respectively. For comparison, the ideal case of ^pk-pk(6') without those effects is 
shown by the dotted lines in the upper left and bottom left figures. The dashed line in the upper left figure 
shows the conventionally used two-point correlation function C{0) normalized to agree with the value of 
Cpk-pk(6') with the solid hne at 6* = 200'. 



structures at scales oi 9 < 9^ as actually shown by the simulated maps. The beam smearing on the intrinsic 
Cpk-pk(^) then appears as a cutoff at scales below the beam size although the effect moderately changes the 
global shape of ^pk-pk(^) in the MAP case (top right and bottom right panels). This is because the beam 
smearing causes an incorporation of intrinsic hotspots contained within one beam. To clarify the noise effect 
explicitly, we also show Cpk-pk(^) only with beam smearing effect (dot-dashed lines). The curves explain 
that spurious hotspots due to the detector noises generate an extra power of correlations on Cpk-pk(^)- In 
particular, the predicted shape of ^pk-pk(^) for MAP is largely affected by the noise effect up to 6* = 80' 
while the Planck cases have slight changes only at 9 ^ 20'. By using equation (A12), we can predict values 
of mean number density of hotspots above = 1 for those cases. The values with both the beam smearing 
and noise effects and only with beam effect in SCDM model are fipk = 8.17 x 10'^ and 7.74 x 10'^[rad"^] 
for Planck (upper left), respectively, while fipk = 5.31 x 10^ and 3.01 x 10^[rad~^] for MAP (upper right), 
respectively. Similarly, in LCDM model ripk = 7.73 x 10^ and 7.33 x 10^[rad"^] for Planck (bottom left) 
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Fig. 5. — The logarithmic power spectra of the gradient (left) and second-derivative (right) fields of the 
CMB anisotropics field as a function of l/l{(x 9). The solid and dashed lines show the results of SCDM and 
LCDM models, respectively. This figure clearly explains that the oscillations of l^Ci is strongly enhanced. 

and npk — 5.06 x 10^ and 2.88 x 10'^[rad~'^] for MAP (bottom right), respectively. These values actually 
coincide with results of the number counts of hotspots in the simulated CMB maps within the Poisson error. 
Moreover, to reveal the differences of shapes between ^pk-pk(^) and the conventionally used two-point 
correlation function of the temperature fluctuations field itself defined by C{6) = {A{0i)A{62))^g^_g^^^g, 
the dashed line in upper left panel shows C{9) normalized to agree with the value of solid line at 6* = 200'. 
It is clear that Cpk-pk(^') has much more oscillatory features than C{6) does (HS; TKF). This reason is as 
follows. In the peak statistics we need statistical properties of the gradient and second derivative fields of 
the temperature fluctuations fleld in order to identify the local maxima (or minima) in the CMB maps. The 
power spectra of those fields per logarithmic interval in I are then I^Ci/{2'k) and l^Ci/ {2n), respectively, 
while the power spectrum of the temperature fluctuations flled is PCi /{2t:) whose integral over I space 
produces C[6). Figure || shows the power spectra of A ^ and A as a function of l/l (ex 6*), and clearly 
explains that they strongly enhance the oscillations of PCi. If recalling the relation of k, 2tt/1, there are 
indeed correspondences between both oscillations of i^pk-pk(^) and CiP/{2tt); the depression a,t 9 ~ 70' 
corresponding to the scales around the flrst Doppler peak, a prominent peak at 6* w 13', and a damping tail 
at 9 < 10' associated with the Silk damping in Ci. Most importantly, these oscillatory features of ^pk-pk(0) 
physically mean that pairs of hotspots are discretely distributed with some characteristic separation angles 
on the last scattering surface. The peak statistics thus produces an additional information that could be 
convenient for the study of the weak lensing although it relys on only a subset of available information (the 
location of peaks, peak height and their profile). Furthermore, we expect that the distribution of hotspots 
are more robust against the systematic observational errors. The previous works already concluded that the 
lensing effect on Ci is small at I < 3000 (e.g. see Seljak 1996 and references therein), and we therefore expect 
that the observed Ci will provide the accurate prediction of unlensed ^pk-pk(^)- The lensing signatures 
to ^pk-pk(^^) are then extracted as non-Gaussian signatures, which are differences between the observed 
(lensed) and the predicted two-point correlation functions. 
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Fig. 6. — The lensed and unlcnsed two-point correlation functions of hotspots as a function of 6 using the 
same parameters of experimental effects as those in Figure ^. The weak lensing effects are computed using 
the same models of matter power spectrum with sets of as as in Figure fl. 



4.3. Lensed two-point correlation function of hotspots 

We present the theoretical predictions of the lensed two-point correlation function of hotspots. By 
using equation we then compute the lensing effect on the intrinsic ^pk-pk(^) into which the experimental 
effects are already included by the method presented in the previous subsection. Although the order of 
computing those secondary effects is not correct, the treatment makes the calculations much simpler and is 
a good approximation because of the following reasons. As explained, the instrumental effects on ^pk-pk(^) 
can be approximately included only by modifying the angular power spectrum C; based on the Gaussian 
theory and actually the predictions are in good agreement with the numerical experiments of ^pk-pk(^) (HS 
and also see Figure This means that the distribution of hotspots as a point process can be accurately 
predicted even in the CMB maps altered by the beam smearing and the detector noise. Then, since we are 
here interested in a problem how the weak lensing redistributes the 'key' hotspots, which can survive after 
the beam smearing effect, in an actual observed map, we can approximately consider the lensing effect on 
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the ^pk-pk(^^) after the instrumental effects as long as the weak lensing does not create a lot of spurious 
hotspots on the map. Moreover, we can at least say that the important lensing signature to ^pk-pk on large 
angular scales such as 6* sa 70' is not directly coupled to the beam smearing effect of Sfwhm ~ 10'- The 
validity of our treatment has been confirmed by the numerical experiments on the relevant angular scales. 
Figure ^ shows both unlensed and lensed two-point correlation functions of hotspots, say ^^k^pk(^) ^^'^ 
Cpk-pk(^), where the threshold v — 1 is similarly assumed and we employ the same models of matter power 
spectrum as in Figure 0. Although = 1 is assumed throughout this paper for simplicity, the two-point 
correlation function of hotspots above height of an arbitrary threshold can be predicted under the Gaussian 
theory and we could use this freedom for measuring the lensing effect as will be discussed later. The upper 
left and bottom left panels for Planck case clearly demonstrate that the weak lensing effect fairly smooths 
out the shape of unlensed ^pk-pk(^) (solid line). The magnitude of the lensing effect is strongly sensitive 
to the amplitude of as in each cosmological models and larger at angular scales where ^pk-pk(^) has more 
oscillatory features (TKF). In particular, we stress that the weak lensing causes a distinct smoothing effect 
on the depression feature of ^pk-pk(^) at scales around 9 « 70'. We therefore expect that the observed 
depth of depression relative to the plateau shape at larger or smaller scales than the scale can be a robust 
indicator of the weak lensing signatures and depends only on the magnitude of ctgl if it is large sufficiently 
to detect. However, since in the MAP cases (upper right and bottom right panels) the detector noise effect 
decreases largely the depth of the intrinsic depression as shown in Figure ^ the lensing effect is hidden 
at the scale for all models. If using the relation between k and / shown in Figure |(a), the measure of 
<7GL {&) from the lensing signature at « 70' can provide a constraint on the amplitude of mass fluctuations 
with large wavelength modes such as A « 50/i~^Mpc. On the other hand, the measurement of ctgl at the 
prominent peak scale such as « 20' is sensitive to the smaller scale structures such as A « 10/i~^Mpc, 
while the shape of fpk-pk is also sensitive to the experimental effects at such scales. 



5. Detectability of the lensing signatures 

In this section, by using numerical experiments we quantitatively investigate how accurately non- 
Gaussian signatures on ^pk-pk(^') induced by the weak lensing can be detected with the expected future 
MAP and Planck surveys. 



5.1. Numerical experiments 

We perform simulations of the CMB maps with and without the lensing effect by using the following 
procedures. A realistic unlensed temperature maps on a fixed square grid can be generated from a given 
power spectrum, C;, based on the Gaussian assumption. Each simulated map is initially on 40 x 40 square 
degree area with 4096 x 4096 pixels. After regridding to take into account the beam smearing effect, the 
actual pixel number is reduced to 1200 x 1200 and 480 x 480 for the Planck and MAP cases, respectively. 
To compute the lensing effect on the CMB maps, we first need the convergence field on each grid. We then 



employ the following angular power spectrum of convergence field ( Blandford et al. 1991 ; Miralda-Escude 
199l| ; [Kaiser 1992|) ; 
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Fig. 7. — An example of simulations of primordial temperature fluctuations map (left) and the map deformed 
by the lensing effect (right) on a side of 2 degree, where we assumed the instrumental resolutions of Planck 
(^fwhm — 5.5' and fiscns = 4.3 x 10~^). We employed the SCDM model and the lensing displacements field 
is computed by using as — 1-5. The peaks above la are marked with crosses and the contours are stepped 
in units of 3asons- 



where the convergence field k is expressed in terms of the radial integral of three dimensional density 
fluctuations field S a.s k — (3/2)HqQ^q JdxW{x,Xrcc)^{x)0'~^^- Note that the second moment of 
convergence field is then ^k^(0))) = J {ldl/2TT)Pf^{l). Using the power spectrum Pk{1), the convergence 
field can be generated as a realization of a Gaussian process. To compute the displacement vector S6^ 
we transform the convergence field in the Fourier space, and compute the Fourier component of the 
displacement vector, 69(1), by using the relation of 

se{i)^2t^kii), (11) 

where we have used the relation of 2k{0) = d{59i) / dOi + d{662) / 862 derived by the cosmological lens 
equation. If transforming it back to real space, for each point on the observed temperature map we 
can obtain the corresponding displacement vector to map the point on a irregular grid of the primary 
temperature map on the last scattering surface. The lensed temperature map can be then obtained by using 
cloud-in-cell interpolation to compute the value on the original regular grid of observed map. In the case of 
taking into account the instrumental effects of beam smearing and noise, we furthermore smooth out the 
temperature map by convolving the Gaussian window function !F{9,6') = 1/ (2^9^) exp[—\6 — 6' \'^ / {29'^)] 
and then add randomly the noise field into each pixel. Figure ^ illustrates an example of realizations of the 
simulated unlensed and lensed maps for the SCDM model with erg = 1.5, which provides the largest lensing 
signatures of our considered models. Interestingly, the figure shows that the displacement of each position 
of hotspots in the lensed map from the unlensed map is relatively large even though the global patterns of 
temperature fluctuations in both maps are not considerably changed. This is the consequence of the large 
scale modes of the deflection angle, which play an important role to the lensing effect on Cpk-pk(^')- 
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Fig. 8. — An example of averaged two-point correlation function with lensing effect expected from Planck 
survey in SCDM (left) and LCDM (right) models with erg = 1.5, respectively. The solid and dashed lines show 
the theoretical predictions of unlcnscd and lensed ■fpk-pk(^) as in Figure ^ and ^, respectively. The nmncrical 
results are obtained from one set of 17 realizations of the lensed maps, and the error bars correspond to the 
cosmic variance computed from 100 realizations. The bottom panels show the results around the depression 
angular scale (w 70') of the unlensed ^pk-pk(^)- 



5.2. The signal to noise ratio of the lensing signatures 

The observational errors associated with measurements of fpk-pk(^) arise from the cosmic variance and 
the instrumental resolutions. To accurately compute the cosmic variance, we have used 100 independent 
realizations of both unlensed and lensed temperature maps, respectively, with 40 x 40 square degree area. 
In this paper we assume 65% sky coverage for MAP and Planck surveys, and this corresponds to the 
assumption that we can use about 17 independent simulated maps for those surveys in order to obtain 
the averaged ^pk-pk(^)- In Figure ^ we show an example of the averaged lensed two-point correlation 
functions of hotspots computed from one set of 17 realizations for the expected Planck runs, where the 
error-bars in each bin can be estimated by rescaling the variance of the estimates obtained from the 100 
realizations by a factor of vTf and the resolution of bins in 6 is limited by the pixel size of simulated CMB 
maps. The unlensed cases have been already presented by HS. The bottom panels also show the results 
around the depression scale 9 « 70' of ^pk-pk(^^) in each cosmological model. Figure ^ clearly shows that 
the theoretical predictions are in remarkable good agreement with the numerical results. In particular. 
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the lensing signatures at the depression scale definitely deviate from the unlensed cases, and we therefore 
expect that the non-Gaussian signatures should be detected by Planck with high significance for some 
adequate values of erg. On the other hand, for the lensing signatures at smaller scales such as the prominent 
scale 9 « 20', it seems to be slightly difficult to distinguish them because the sampling variance are larger 
at smaller scales and the shapes of Cpk-pk(^) is also sensitive to the instrumental effects of beam smearing 
and noise as shown in Figure ^. On the other hand, we concluded that it is difficult for MAP survey to 
distinguish the lensing signatures on ^pk-pk(^) for all cosmological models considered in this paper. For 
this reason, we present only results from Planck in the following. 

We so far have adopted the specific shape of matter power spectrum in each considered cosmological 
model, and therefore a free parameter of our model is only as- While the actual observable quantity of 
our method is the lensing dispersion ctgl (9) , we here investigate the dependence of the lensing signatures 
to Cpk-pk(^) on (Tg for simplicity and the problem will be discussed in the final section. In Table ^ we 
summarize the results obtained for the determination of ug with a best fit and the error associated with this 
determination. The error then has been determined as the variance of the best fit values of erg among a lot 
of sets of the numerical experiments performed by fitting the simulated results to the theoretical templates 
of lensed ^pk-pkl^*) so that the x-square is minimum. Each best fit of theoretical curves to the numerical 
experiments is restricted mainly from the data at > 20', in particular the data around the depression scale 
of ^pk-pk(^)- One can see that the signal to noise ratio indeed grows with erg in each cosmological model. 
Furthermore, if comparing results obtained from SCDM and LCDM models with the same value of cg 
(ug ~ 1.5), it is clear that more robust constraints can be provided in SCDM, more generally for background 
models with larger Q^no, as expected. The noise level of Planck is independent on the erg estimations, but 
we have confirmed that the beam size is rather important for the detections. Importantly, the accuracies 
of (Tg determinations expected from Planck reach about 8% and 11% for models with cg — 1.5 in SCDM 
and LCDM, respectively. Furthermore, in the Gaussian random theory the two-point correlation function 
of local minima [coldspots) should have the same shape as that of hotspots, and therefore combing the 
measurements of coldspots correlation function improves the lensing signals by about factor We could 
also improve the significance by combining the measurements of ^pk-pk with another different thresholds, 
but the independence of data then has to be carefully investigated because the angular positions of same 
hotspots are used many times in the fitting. 

Our arguments presented in this section rely on the expectation that we can accurately predict the 
shape of unlensed ^pk-pk(^) as a function of sets of cosmological parameters constrained from the measured 

Table 1: Summary of best fit of txg determinations from the lensing signatures to ^pk-pk(^) using the 
numerical experiments for the expected Planck survey in SCDM and LCDM models. We have assumed 65% 
sky coverage, and the Icr error in each determination represents uncertanities caused by the cosmic variance 
and the instrumental effects due to the beam smearing effect and the detector noise (see text). 



input values of erg 


SCDM (a„o = l,/i = 0.5) 


LCDM (r^mo = 0.3, Vlxo = 0.7, h = 0.7) 





0.30 ±0.31 


0.23 ±0.35 


0.5 


0.47 ±0.24 




1.0 


1.03 ±0.13 


0.93 ±0.23 


1.5 


1.52 ±0.12 


1.49 ± 0.17 


1.5 (without noise) 


1.47 ±0.12 


1.53 ±0.16 


2.0 




1.97±0.16 
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Fig. 9. — This figure shows the fake theoretical predction of Icnsed hotspots correlation function (solid line), 
^pk-pk, computed from the lensed (measured) Ci based on the Gaussian theory in the SCDM model with 
(Js — The exact prediction (dashed line) and the unlnesed case (dotted line) are also shown. Since 
the lensed temprature fluctuations field weakly deviates from the Gaussian, ^pk-pk can no longer accurately 
describe the distribution of hotspots in the lensed maps. 



angular power spectrum Ci within the limit of the observational errors. However, we should bear in mind 
the fact that only the lensed C; is measurable. Then, one may imagine an approach to compare the 
measured ^pk-pk to the fake prediction, say ^pk-pk, computed directly from the lensed (measured) Ci 
based on the Gaussian theory. Since the lensed temperature fluctuations field weakly deviates from the 
Gaussian (Bernardeau 1997), the distribution of hotspots in the lensed maps can be no longer characterized 
only by the lensed C;. However, it will be still interesting to see differences between the exact and fake 
predictions of the lensing effect on ^pk-pk(^)- Figure ^ thus shows the shape of ^pk-pk and reveals that 
Cpk-pk overestimates the power of correlations on scales oi 9 < 40', while ^pk-pk mimics the smoothing 
effect on the depression feature of ^Jpk-pk resulting from the fact that the weak Insing already causes the 
smoothing of C;. Consequently, for example, the value of reduced x-square for the fitting between ^pk-pk 
and the simulated ^pk-pk becomes worse to be ~ 162/95 against the fact that the fitting of using 
the exact theoretical predictions @j of lensed ^pk-pk produces values around almost unity. Hence, the 
important problem how we extract the lensing information only from the measurable CMB quantities has 
to be carefully investigated and this problem will be again discussed in the next section. 
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6. Discussion and Conclusions 

In this paper, we have quantitatively shown that the peak-peak correlation function for the CMB maps 
can be an efficient measure to probe the mass fluctuations at relatively large scales. The non-Gaussian 
signatures on ^pk-pk(^) caused by the weak lensing appear at several scales as smoothing effects on the 
oscillatory features of ^pk-pk(^) (TKF). In particular, by using numerical experiments of the CMB maps 
including effects of beam smoothing and detector noise as well as the lensing displacements, we found 
that the lensing signatures at about one degree scales, where the intrinsic ^pk-pk(^) has the pronounced 
depression feature, should be detected most significantly from the expected Planck survey if the signatures 
are adequately large for the detections. On the other hand, unfortunately we concluded that it is difficult 
for MAP to detect the lensing signatures for the cosmological models considered in this paper. The direct 
observable quantity of our method is the dispersion of lensing deflection angle as shown by equation (^. 
We have then revealed that the one-degree scale dispersions are sensitive to amplitudes of three dimensional 
mass fluctuations around wavelength A ~ 50/i~^Mpc because of the projection effect. Furthermore, such 
lensing dispersions have contributions from structures of the universe with a wide redshift distribution in 
the range of 1 ;S z ;S 3 (see Figure This projection effect would be thus a serious issue for extracting 
the cosmological information from weak lensing contributions in a general case. However, since the mass 
fluctuations around A ~ 50/i~^Mpc are now still in the linear regime and the evolution history in the 
redshift space is theoretically well understood in the context of gravitational instability in an expanding 
universe, the lensing signatures can accurately determine the cosmological parameters associated with 
amplitude and evolution of mass fluctuations at the scales. It is therefore expected that our method can 
provide robust constraints on cosmological parameters cts a-nd flmo without much specifying the shape of 
matter power spectrum. Our numerical experiments indeed revealed that significant signal to noise rations 
for determinations of as from the lensing signatures to ^pk-pk(^) are obtained for some input values of erg 
in SCDM and LCDM models, respectively. It was also shown that for the same value of as more significant 
signal can be obtained in SCDM than in LCDM. To break the degeneracy in as — ^mo determinations 
only by using our method, one must measure the scale dependence of lensing signatures. This seems to be 
difficult even for Planck survey, because the shape of ^pk-pk(^) at small scales of ^ ~ 20' where ^pk-pk has 
the secondly significant lensing signatures is also sensitive to the instrumental effects of beam smoothing 
and noise. Anyway, it is very interesting that our method could provide constraints on the mass fluctuations 
in the linear regime independently of those provided by the survey of galaxies clustering and the primary 
CMB anisotropies alone, which cannot directly probe the power spectrum of dark matter. 

In the results presented in this paper, we have discussed how accurately the non-Gaussian signatures 
of the lensed C^i^^lpj^ can be detected as a deviation from the unlensed ^pk-pk as shown in Figure ^. This 
strategy is based on the assumption that we can accurately predict the unlensed fpk-pk(^) from the 



measured C; because the lensing effect on Ci is small ( Seljak 1996b ). However, since an actual measurable 
quantity is only the lensed C;, there remains an important problem that we have to carefully investigate. 
Usually, the measured angular power spectrum C; can be used to constrain the sets of cosmological 
parameters under a specific scenario, for example, within the framework of inflationary-motivated models 
(e.g. Lange et al. 2000| ). The detailed analysis of using higher / modes such as / > 3000 will need to 



take into account the lensing contributions to Ci for the accurate determinations (Stompor & Efstathiou 



1999| ). Therefore, one approach toward detecting the lensing effect on fpk-pk is that we first construct a 



lot of templates of the unlensed ■^pk-pk as a function of sets of cosmological parameters constrained by the 
measured Ci within the observational errors, and then we compare the templates with the directly measured 
Spk-pk in the observed (lensed) CMB maps by taking into account the lensing contributions. Based on 
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the considerations, we are now investigating the detailed dependence of cosmological parameters on the 
shape of ^pk-pk- As a resuh, we confirm to a extent that, if we focus on the depth of depression feature of 
Cpk-pk(^) at one degree scales relative to the plateau shapes at larger or smaller scales than that scale, the 
measured depth could be a robust indicator relatively independently of the cosmological parameters because 
the weak lensing with same magnitude of ugl can shallow the measured depth to a same amount. This 
should be further investigated carefully and will be presented elsewhere. The important thing that we stress 
in this paper is that we quantitatively proposed a new statistical method for measuring the lensing effect 
based on the peak statistics. It is not evident that our method and the other methods can measure the 
small lensing signatures from the observed CMB data with exactly same statistical significance. Therefore, 
several independent methods should be performed to measure the lensing signatures complementarily. 

Undoubtedly, we have to carefully consider how secondary anisotropics and the foreground sources affect 
our conclusions. The most important source of secondary anisotropics is the (thermal) Sunyaev-Zel'dovich 
(SZ) effect, which is essentially caused by hot electrons inside the clusters of galaxies. Since SZ induces 
redistributions of the photon energy from low frequency to high frequency part of the black body spectrum, 
the effect could mimic peaks in the observed temperature map. However, the SZ anisotropics are dominated 



by the Poisson contribution from the individual clusters at relevant angular scales (Komatsu & Kitayama 



1999| ), and therefore the contribution of the cross correlation between spurious peaks due to SZ and the 



intrinsic peaks will be smaller than the amplitude of the intrinsic peak-peak correlation because the peaks 
on the last scattering surface and the SZ clusters are statistically uncorrelated. Furthermore, we should 
emphasize that the SZ effect can be always removed by either observing at 217 GHz or by taking advantage 
of its specific spectral properties. The other secondary anisotropics such as Rees-Sciama effect do not affect 
our method because the amplitudes of those effects have a very small contribution at I < 3000 ( peljak 



1996a). Finally, with regard to the effect of the extragalactic sources, since it is expected that the amplitude 
of anisotropics due to the discrete sources in the 100 — 200 GHz range are well bellow the amplitude of 
primordial fluctuations ( Toffolatti et al. 199§| ) and the sources will be also eliminated by using the multi 



frequency observations in principle, it can be safely concluded that this is not a serious for our method. 



Recent measurements of cosmic shear (e.g. see Van Waerbeke et al. 2000| ) have provided constraints 



on the combination of tig and flmo- The cosmic shear can probe the mass fluctuations at angular scales 
of 1' ;S 6* ^ 30' because of the limited survey volume, where the nonlinear clustering effect of dark matter 
could play an important role. It has been also shown that the deformation effects on the CMB maps could 



be measured by Planck ( |Bcrnardeau 1998| ; |Waerbeke, Bernardeau fc Benabed 1999| ), and it is sensitive 



to the projected mass fluctuations at scales around 5' which is the characteristic curvature scale of the 
primary temperature field. Therefore, although flmo parameter cannot be determined accurately form the 
primary CMB anisotropics alone even with Planck because of the so-called cosmic (geometrical) degeneracy 



(Bond, Efstathiou, & Tegmark 1997), it is expected that our method and those other independent methods 
of using the weak lensing will break the cosmic degeneracy with high precision because the amplitude of 
lensing contributions is also sensitive to ilmo as explained. Another challenging possibility is that those 
independent methods could allow us to observationally reconstruct the shape of power spectrum of dark 
matter including the evolution history in the redshift space by combining those measurements of amplitudes 
of mass fluctuations at respective angular scales. 
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A. Mean density of hotspots in the 2D CMB map 

In this appendix, we briefly review the derivation of the mean number density of peaks in the two 
dimensional Gaussian field ( Bardeen et al. 1986 , Bond fc Efstathiou 1987 ). 

Following BBKS and BE, we introduce the spectral parameters defined by 



7 ^ 



V2^, (Al) 



ao(T2 (72 

where 

- f-icd'-. (A2) 

The characteristic curvature scale of primary temperature field can be estimated as 9^, ~ 5' for the 
cosmological models considered in this paper. 

The problem of the two dimensional peak statistics is to consider the statistical properties of the 
point process. We can then express the point process entirely in terms of the field and its derivatives. 
At hotspots (coldspots) the gradient Ai{= d/S./dO'^) vanishes, and the eigenvalues of curvature matrix 
Aij(= A/ {89^86^)) are all negative (positive). We therefore need to consider six independent variables 

V = {A^ Ay, Axx, Ayy, Axy) to specify one local maximum. For the Gaussian field, the probability 
density function (PDF) for v is 

= (2.)3|detW)l^/^ ' (A3) 

where the covariance matrix is Af^ — (viVj) because of (vi) = in the present case, and is the inverse 
matrix of Mij. Following HS, we introduce the notations for convenience as 

_ A _ A 

v= — , 7]i = — (i = a;,y) 

(To <Tl 

_ AxX + Ayy ^ _ AxX ~ Ayy ^ _ 2Axy (A4) 

(72 (T2 (T2 

The non-zero second moments of these variables are 

(i.2) = (A2)=2(r2)=2(Z2)=2(r;f) = l, (^A) = 7, (A5) 

where 7 = (7i/((To(72)- By using these equations, we can obtain the following PDF for the variables 

V = (v, A, Y, Z, rji) with the simple form, and it gives a probability that the field point has values in the 
ranges oi v to v + dv, A to A to A + dX and so on: 

92 

Pi{v,X,Y,Z,'ni)dvdXdYdZd^ri, = Pvp [~Q] di,d XdYdZd^ri, (A6) 

(27r)3^1 - 72 
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with 



2Q = 



(A7) 



1-72 

The point process of hotspots (or coldspots) can be described by the number density 'operator' 

"pkW =5I'5''(^-^pk.p), (A8) 



where 5^{x) is the Dirac delta function. In the neighborhood of a hotspot point O^y we can expand the 
field A(0) in a Taylor series: 



(A9) 



Using this equation, therefore, the number density field can be expressed in terms of the Gaussian variables: 

7ipk(0) - |det(A„(0))|(5^(A,(0)). (AlO) 



The summation symbol in equation (A8) can be eliminated because the delta function of the above equation 
picks out all of the extremal points which are zero of Ai(0) are maximum in the two dimensional map. 



Hence the ensemble average of equation (AlO) produces the differential mean number density of 
hotspots of height in the range of v and v + dv. 



1 



n^^iu) = {\detiA,,)\5^{A,)) = ^{\X^-Y'-Z^\6^iri,)) 



1 

2^ Jo 



dX I dY 



dZ{X^ -Y^ - Z^)pi{u, X, Y, Z, 77, - 0), 



(All) 



where we have adopted the conditions for a hotspot of X > and Y^ + Z^ < X^. BE derived the analytical 
expression for f;,pk(j^); 



1 



(27r)3/202 



exp(-i.V2)G(7,7f^), 



(A12) 



where 



G(7,x,) = ixl^^^){l--CTic 



+ 



exp[-x2/(3-272)] 
(3-272)1/2 



v/2(l-7^) 



1 erfc 

2 



+ a;,(l-7-') 



2,exp{-x2/[2(l-7^)]} 



(A13) 



y2(l-72)(3- 272) 

In this paper, we have often used the mean number density of peak of height above a certain threshold 
'^thresh obtained by integrating equation ( A12 ) over v > I'thrcsh- Naturally, the mean number density of 
coldspots of height below —v (Apk < —i^ctq) is symmetrically given by ncoid(< ~i^) ~ ^pk(> i^)- 



B. The unlensed two-point correlation function of hotspots 



Based on the peak statistics for the two dimensional Gaussian field, we briefly review the derivation 
of two-point correlation function of hotspots following HS. For this purpose, let us consider two hotspots 
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separated by the angular scale Q. In the same way as in appendix Q we then need the joint probability 
density function for the 12 independent variables v = (vx^v-i) = (y(Q\)^v(d2)) with Q — \B\ — 02|- If using 
the variables (A4) for each hotspot, we can block the covariance matrix for the 12 variables in the order 
(i^i, Xi, ryia;, 5^1, 1^2, -^^2) ??2£!;j '7ii/j f]2y^ Zi, Z2) into the 8 X 8-matrix Mjg-) (upper left) and the 4 x 4-matrix 
M(4) (bottom right); 
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M(4) = 
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where we have introduced the notations (HS): 
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Ztt 



(Bl) 



(B2) 



(B3) 



Note that both Mg and M4 are symmetric. 

The joint PDF for the two hotspots can thus be expressed in terms of the above covariant matrix; 



1 



(27r)6|det(A%)P/' 



exp 



1 



(B4) 



Therefore, as in the similar way of the derivation of mean number density, we can calculate the unlensed 
two-point correlation function of hotspots which are separated by 9 and have heights above vi and V2, 
respectively; 



1 + Cpk-pk(^) 



?T.pk(> J^l)?T.pk(> V2) 
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46'Jnpk(> vi)np\^{> rt2) 
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pk(0l)'^pk(^2))|0^„0^| = 

((xf-yf -zf)(x 
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^2 



dX'2 



z'i)s^{v'u)s''{v'2.)),e._e^^=e 



dY[ 

-XI 



dYi 



-X' 



dz[ 



dZ'2{X[^-YC-ZT)(X'i-Yi 



^P2{v'l,X[,Y{,Z[,J^[^ = 0,V2,X2.Y2,Z2,'q'2i = 0). 



(B5) 
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To obtain ^pk-pk(^), we performed a six dimensional numerical integration of the equation ( [B5| ) because 
the two integrals can be done analytically. 



C. Alternative representation of the lensed hotspots correlation function 



The lensing dispersion ctgl(^) generally has isotropic and anisotropic contributions, which are expressed 
by o'GL,o(^') and (Tgl,2(^), respectively, as shown in equation (||). For the standard power spectra of mass 
fluctuation as adopted in this paper, the contribution of (Tgl,2(^) is smaller than that of crGL.o(^)- If the 
neglecting the contribution, the key equation of this paper can be further simplified. In this case, the 
equation (ph becomes 



de'e' //d/epk-pk(e')^o(^f?')exp 



Cpk-pk(6') R 

and this can be then rewritten as 

epk-pk(e) - l0'd9'Cp^.-pi.{0')K{e,9'), 
where the kernel K{9, 6') is given by 



K{e,e') 



ldlJo{W)Jo{W')exp 



f^GL,o(^) 



exp 



Q2 ^ Q>2 - 




ee' 










.^^L,o(^). 



(CI) 



(C2) 



(C3) 



/o(a;) is the modified zeroth-order Bessel function, and we have used equation (6.633.2) of Gradshteyn & 
Ryzhik (1994). As shown in Figure |l|, (TQi^fi{e)/e < 1 and the argument of Iq is generally a large number. 
Therefore, if using the approximation Io{x) ~ (27ra;)^^/^ exp(x) for x — > oo, we can rewrite equation ( |Cl| ) 
in the following form; 



GL 

pk— pk 



(0) 



{2^eY/^aGi.,m 



e'de'^pk-pk{e') exp 



Qf\2 



(C4) 



This expression explicitly explains that the lensing effect on ^pk-pk(^) appears as a Gaussian smoothing with 
relative width (Tgl.Oj and then the asymptotic behavior at (Jgl,o — > is naturally C^k^pk(^) Cpk-pk(f^)- 
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